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^ We study the collider phenomenology ol bilinear R-parity violating supergravity, the simplest 

effective model for supersymmetric neutrino masses accounting for the current neutrino oscillation 
data. At the CERN Large Hadron Collider the center-of-mass energy will be high enough to probe 
' directly these models through the search for the superpartners of the Standard Model (SM) particles. 

' We analyze the impact of R-parity violation on the canonical supersymmetry searches - that is, we 

examine how the decay of the lightest supersymmetric particle (LSP) via bilinear R-parity violating 
interactions degrades the average expected missing momentum of the reactions and show how this 
diminishes the reach in the usual channels for supersymmetry searches. However, the R-parity 
violating interactions lead to an enhancement of the final states containing isolated same-sign di- 
leptons and trileptons, compensating the reach loss in the fully inclusive channel. We show how 
the searches for displaced vertices associated to LSP decay substantially increase the coverage in 
supergravity parameter space, giving the corresponding reaches for two reference luminosities of 10 
and 100 fb~^ and compare with those of the R-parity conserving minimal supergravity model. 
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I. INTRODUCTION 



Weak scale supersymmetry (SUSY) represents the most popular approach towards a solution of the 
hierarchy problem, with the added virtue of being perturbative and accounting for the unification of 
the gauge coupling constants and the possibility of eventually including gravity Little is currently 
known about the detailed mechanism of how SUSY is realized and of how it is broken. Ultimately 
only experiment can settle this issue. The most popular ansatz assumes that R-parity is conserved 
and that SUSY breaking in a hidden sector is communicated to the observable sector by flavour-blind 
gravitational interactions. Neither of these assumptions is mandatory, in fact they are both ad hoc. Many 
different mechanisms of supersymmetry breaking can be envisag ed 0, S, 3 and, similarly, the breaking 
of R-parity Q can take place in a variety of ways @, 0, S] ■ 

Here we will adopt only the first part of the conventional ansatz, namely assume transmission of super- 
symmetry breaking via gravity but relax the R-parity conservation assumption. A specially interesting 
possibility is that R-parity breaks spontaneously, as a result of the minimization of the scalar Higgs po- 
tential [9|, very much in the same way as the electroweak symmetry itself breaks. Spontaneous R-parity 
breaking models are characterized by two types of sneutrino vacuum expectation values (vev) 48| , those 
of right and left-handed sneutrinos, singlets and doublets under SU(3) ^ SU(2) (E) U(l) respectively, 
obeying the "vev-seesaw" relation vlVr ~ h^rri^ where hi, is the small Yukawa coupling that governs 
the strength of R-parity violation [9] . Below the scale vr the effective breaking of R-parity is explicit 
through bilinear superpotential terms and the corresponding soft supersymmetry breaking terms. For 
generality and simplicity here we replace the right-handed sneutrino vr vev by the effective bilinear 
coupling. 

We refer to such simple effective bilinear R-parity breaking (BRpV) picture as BRpV-mSUGRA for 
short [l3| [ill ■ Here we start directly from the BRpV-mSUGRA model, which has a strong phenomeno- 
logical motivation, besides simplicity, namely that it provides the most economical effective model for 
supersymmetric neutrino masses, that successfully accounts for the observed pattern of neutrino masses 
and mixings [l^, E3|- Moreover, it has been shown that, despite the smallness of R-parity violation 
indicated by the neutrino oscillation data, its effects can be tested and the model falsified at collider 
experiments [l3| ■ Here we focus on the implications of the model for the first stages of the CERN Large 
Hadron CoUider (LHC). 

There have been extensive studies of searches for supersymmetric particles but mainly in the context of 
models with conserved R parity 
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16, 17 



So far no direct experimental evidence for SUSY in high 
energy experiments has been found. With the coming into operation of the LHC searches for R-parity 
"3, 21, 22, 23 1 acquire stronger motivation not only by their intrinsic interest. 
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breaking signatures 

but also because they o pen the exciting possibility of probing the neutrino oscillation angles at collider 
experiments [2J, bsj, l2a bj, ba, bd, [sO | . 
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Previously we have studied the Tevatron potential to search for BRpV-mSUGRA in the multilepton 



channel [l^ as well as in the form of displaced vertices 20 1. In this work we provide a quantitative 
analysis of the phenomenological implications of bilinear R-parity breaking supergravity for the LHC, 
extending our studies to the new energy range accessible at the LHC and including all constraints from 
neutrino oscillation physics. Since the effective strength of R-parity violation must be small in order to 
fit current neutrino oscillation data, our study also constitutes a robustness check of the supergravity 
parameter reach estimates against the presence of "perturbative" BRpV terms. Supersymmetric particle 
spectra and production cross sections are expected to be the same as the conventional ones, thanks to the 
smallness of R-parity violation required by neutrino oscillation data. Similarly, processes such as 6 — > 57 
and g-2 are essentially the same in BRpV-mSUGRA as in mSUGRA and hence the resulting constraints 
for the latter still hold. 

The basic difference between BRpV-mSUGRA and the conventional R-parity conserving mSUGRA 
scenario is that the lightest supersymmetric particle is no longer stable and it decays, typically inside 
the detector, consequently modifying the SUSY phenomenology at colliders. This has two important 
implications: 

• the neutralino is no longer a viable dark matter candidate; one needs to imp lement either an axion- 



like, gravitino 3l|, |33, l33|, |3J| or possibly majoron dark matter scheme 
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361. Clearly no "SUSY 



dark matter" constraint should be applied in these models. 

• any SUSY particle can be the lightest one, hence it can be electrically charged. For the case of 
interest assumed here it means that the parameters that lead to stau-LSP are perfectly viable for 
the searches, and are indeed included in our analysis. 

As far as searches go, the LSP decay implies that the missing transverse momentum of the events gets 
diminished and, correspondingly, an increase of the multiplicities for jets and/or charged leptons in the 
final state is expected. Most remarkably, within a large range of BRpV-mSUGRA parameters, the LSP 
can live long enough to give rise to displaced vertices, opening a new window for its search and study. An 
interesting phenomenological feature of BRpV not present in generic RPV schemes is that it naturally 
leads to displaced vertices 49]. 

In this work we have analyzed in detail the impact of the LSP decay on the potential of the LHC to 
unravel the existence of supersymmetry in the simplest framework where its breaking is gravitational. To 
achieve this we start from the same conventional R-parity conserving search channels. We demonstrate 
that the LHC reach is reduced in the jets and missing transverse momentum channel with respect to 
the R-parity conserving scenario. Furthermore, we also find that the final state topologies containing 
isolated charged same-sign di-leptons and trileptons are enhanced, partially compensating the loss of 
missing momentum brought about by the LSP decay in the all inclusive channel. An especially relevant 
fact is that the LSP decay can also lead to displaced vertices. We have also investigated the capabilities 
of the LHC experiments, ATLAS and CMS, to unravel the existence of SUSY with R-parity violation 
through the search for displaced vertices with a high invariant mass associated to them. In fact, we 
find that at 100 fb^^ the reach in the displaced vertices is typically larger than the one expected in the 
R-parity conserving scenario, especially at large mo values. 



II. THE BRPV-MSUGRA MODEL 



As already mentioned, here we focus on the phenomenology of the effective theory described by the 
bilinear R-parity breaking supergravity model in which the MSSM superpotential is supplemented by the 
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following terms [lot 

WbrpV = Wmssm + eab^iL^Hl^ , (1) 

with three extra parameters {ci), one for each fermion generation. In addition to these bilinear terms, we 
must also include new soft supersymmetry breaking terms {Bi). 

Kioft = ^MSSM ~ £ab [Bi^iLIH^] . (2) 



Taken together, the new terms in the BRpV Lagrangian (the three bilinear and the Bi) lead to the 
explicit violation of Icpton number as well as R-parity. These terms also induce vevs Vi = vn, i — 1,2,3 
for the three left-type sneutrinos. Note, that in general there is no basis where both sets of bilinear 
R-parity breaking terms can be eliminated at the same time. In the basis where the bilinear terms in the 
superpotential are rotated away one would obtain a model where the trilinear terms are approximately 



given by A-yfe = (ei//i)^£;jfc and A ■ = {ei/^i)YD,jk [37 1. 



To further specify the model we assume a minimal supergravity framework with universal soft super- 
symmetry breaking terms at the unification scale. This model contains eleven free parameters, namely 

Too , mi/2 , tan^ , sign(//) , , , and , (3) 

where mi/2 and mg are the common gaugino mass and scalar soft SUSY breaking masses at the unifi- 
cation scale, ^0 is the common trilinear term, and tan/3 is the ratio between the Higgs field vev's. For 
convenience, we trade the soft parameters Bi by the "alignment" parameters — CiVd + /ifj which are 



more directly related to the neutrino-ncutralino properties 13| . 

The bilinear R-parity violating terms give rise to mixing between Standard Model and SUSY particles. 
For example, in these models neutrinos mix with the neutralinos giving rise to the neutrino masses and 
mixings after the diagonalization of a 7x 7 mass matrix for neutrinos and neutralinos. At the tree level only 



one neutrino picks up a mass while the others acquire masses only radiatively |12l . Il3l [IJ, l38l | . Typically 
the tree-level scale is the atmospheric scale and, at this approximation, solar neutrino oscillations do not 
take place. The solar mass splitting is "calculable" and the solar angle acquires a meaning only when 



loops are included. It has been checked that current neutrino oscillation parameters [39| can be well 



reproduced provided |ei| <C where /i denotes the SUSY bilinear mass parameter [13 1. 

The LSP lifetime depends upon the SUSY spectrum and the R-parity violating parameters. In order 
to satisfy the neutrino constraints, the R-parity violating interaction must be feeble, thus leading to a 



lifetime large enough to give rise to displaced vertices [27[. For the sake of illustration, we depict in 
Figure [1] the Xi decay length for 100 GeV < mo < 2000 GeV, 300 GeV < mi/2 < 800 GeV, = -100 
GeV, tan(/?) — 10, and sgn(/i)=-|-l. We scan randomly all possible BRpV parameter values ei and A^, 
and for each SUGRA point in the plot, we accumulate hundreds of solutions compatible with neutrino 
oscillation data at 3(t level 39|. In this way the top (red) band corresponds to toi/2 = 300 GeV, and all 
possible values of Ci and A^ compatible with neutrino mixings and squared mass differences at 3f7. 

It is important to notice that in our BRpV-mSUGRA scheme, for a given mSUGRA parameter there 
is a range of variation in the BRpV parameters and Aj. While this has no effect in the production 
cross sections of supersymmetric states, it affects the LSP decay length, as illustrated in Fig. [1] However, 
the decay length can vary by not more than ~ 30% for a fixed value of toi/2 over the allowed range of 
the BRpV parameters, except at small toq where light scalars play an important part in the LSP decay. 

Apart from generating neutrino masses, neutralino-neutrino mixing also leads to decay of the LSP into 
Standard Model particles. In scenarios where the lightest neutralino is the LSP its main decays are 

• leptonic decays: 
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200 400 600 800 1000 1200 1400 1600 1800 2000 

mo [GeV] 

FIG. 1: xi decay length versus mo for Aq — —100 GeV, tan/3 — 10, fi > 0, and several values of mi/2. The 
widths of the three shaded (colored) bands around mi/2 = 300, 500, 800 GeV correspond to the variation of the 
BRpV parameters in such a way that the neutrino masses and mixing angles fit the required values [s^ within 
So-. 



- Xi -> with £ = e, ^ denoted by U; 



Xi 



called 



semi-leptonic decays: 



X? 
X? 



vqq denoted jj\ 

Tq'q, called rjj; 

iq'q, called £jj; 

vbb, that we denote by bb; 

vbb, that we denote by bb; 



• invisible decays: Xi ^ vvv. 

We should mention that many intermediate states contribute to the above decays, hke gauge bosons 
or scalars. Depending on the spectrum, some of these final states can be dominated by 2-body decays. 
For instance, the decay Xj' vt^t~ can be dominated by the 2-body Xi — > vZ followed by Z ^ . 
Another example is the 2-body mode vh with h bb. Consequently we single out the 2-body 

decays when presenting our results. 

Figure [2] presents separately neutralino branching ratios into two-body and three-body decay final 
states as a function of mo for Aq = —100 GeV, tan/3 = 10 and /i > 0. The top panels correspond to a 
fixed value of toi/2 = 400 GeV, and the bottom one correspond to mi/2 =800 GeV. In the left panels we 
give branching ratio predictions for optimized and Aj BRpV parameter values that reproduce the best 
neutrino oscillation parameters [39 1. 
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FIG. 2: Lightest neutralino branching ratios as a function of mo for = —100 GeV, tan/3 = 10, and > 0. In 
the upper panels we have mi/2 = 400 GeV while the lower ones the correspond to m-1/2 ~ 800 GeV. Explanation 
in text. 



Recent neutrino oscillation data from reactors (KamLAND) and accelerators (MINOS) have brought 
an improved determination of the solar and atmospheric mass squared splittings [39||. On the other 
hand the solar angle is relatively well determined by current solar neutrino data. In contrast, current 
uncertainties in the determination of the atmospheric angle remain large and unlikely to improve in the 
near future. This translates into a large uncertainty in the expectations for the LSP decays into gauge 
bosons, as illustrated in the right panels in Fig.[2l 

The light and dark-shaded grey bands in the right panels in Fig. [5] indicate, as examples, the uncer- 
tainties for Xi ~^ W~^T~ and Xi ~^ vt~^t~ respectively (note that similar uncertainties exist for other 
channels. For example, x5 ~^ ^~ has a similar spread as — > W~^t~ but, to avoid an over-crowded 
plot, we have not displayed). The width of these bands illustrates the effect of varying the BRpV param- 
eters randomly while satisfying the current neutrino oscillation data, especially the atmospheric angle. 



7 



at 3(7. 

As one can see, for best-chosen values (left panels) the main LSP decay mode at small mo is into lepton 
pairs including a tau, due to light scalar tau exchange. In contrast, for large mo values the main decays 
are into charged lepton (e, fi or r) accompanied by jets coming from a W. 



III. ANALYSIS 

In what follows we will describe our two main analysis methods employing first the conventional SUSY 
search methods and then taking advantage of the displaced vertex technique. As already mentioned 
all the SUSY particle decays, both R-parity violating as well as R-parity conserving, are calculated 



using a generalized version of the SPheno code 40| 50|. Given a point in SUGRA parameter space 



SPheno searches for a set of Ci and A^ that leads to neutrino masses and mixings compatible with current 
neutrino oscillation data. In general, there is a range of R-parity violating parameters that satisfy the 
above neutrino constraints. As already mentioned, both SUSY spectra and production cross sections are 
well approximated by the conventional R~parity conserving mSUGRA ones, in view of the smallness of 
R-parity violation parameters required to fit the neutrino oscillation data. However, for consistency, our 
version of SPheno calculates spectra and branching ratios including BRpV effects. 



We employed PYTHIA version 6.409 [4l|, |42| to generate the signal and backgrounds, inputting the 
SPheno output in the SLHA format |43|. Note that we have considered all SUSY production processes 
in PYTHIA since they all contain two LSP's in the decay chain. Moreover, we have assumed that the 
detector possesses cells 0.1 x 0.1 in A(/) x A?/ and included a Gaussian smearing of the energy with 
A(£') — 0.5 X \/E {E in GeV). Jets were defined using the subroutine PYCELL with a cone size of 
AR = 0.7. We chose the CTEQ5L parton distribution functions. 



A. Standard searches 

The LSP decay inside the detector certainly has an important impact in the supersymmetric particle 
searches at the LHC. For example, visible LSP decays reduce the missing transverse energy leading to 
a weakening of the standard supersymmetry signals. In order to access the impact of the LSP decay in 
the usual SUSY searches we analyzed the LHC potential for searching for SUSY in the same final states 



analyzed in [44| , i. e. 



1. Inclusive jets and missing transverse momentum (denoted by (IN)): in this case just jets and 
missing are taken into account in the analysis; 

2. Zero lepton, jets and missing transverse momentum (denoted by (0^)): here we consider only events 
presenting jets and missing without isolated charged leptons (e* , A^^); 

3. One lepton, jets and missing transverse momentum (denoted by {11))'- we include in this final 
state topology only events presenting jets and missing -^t accompanied by just one isolated charged 
lepton (e^ , m"*")- Notice that final state r's contribute to this topology only through their leptonic 
decays; 



4. 



Opposite sign lepton pair, jets and missing transverse momentum (denoted by (OS)): this final state 
is characterized by the presence of jets, missing ^t, and two isolated leptons of opposite charges; 
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5. Same sign lepton pair, jets and missing transverse momentum (denoted by (SS)): here we consider 
only events presenting jets and missing -^t accompanied by two isolated leptons of the same charge; 

6. Trileptons, jets and missing transverse momentum (denoted by (3£)): we classify in this class events 
that exhibit jets, missing and three isolated charged leptons; 

7. Multileptons, jets and missing transverse momentum (denoted by (M^)): this final state presents 
jets and missing accompanied by more than three isolated charged leptons. 



B. Standard event selection 



Our analysis of the topologies defined above follows closely the procedure of 4J]. Initially, we applied 
the acceptance cuts: 

ACl We required at least tvifo jets in the event with 

> 50 GeV and \rij\ < 3 . (4) 

AC2 In order to suppress a large fraction of the enormous two jet background, we demanded that the 
transverse sphericity of the event exceeds 0.2 

St > 0.2 . (5) 

ACS To reduce the background stemming from mismeasured jets, the missing transverse momentum 
should not be aligned with any jet. Therefore, we imposed that the azimuthal angle (Ayj) between the 
jets and the missing momentum must comply with 

g<A^<-. (6) 
In the studies of the different final state topologies we used a floating cut designed to give some 



optimization of the cuts over the parameter space 4J|. We considered = 200, 300, 400, and 500 GeV 



Given one value of E^, we required for all final states that 



p^'^ > E^ and i)T > E^ (7) 

where pl^"^ stand for the transverse momenta of the two hardest jets. These are the only cuts applied 
for the IN topology in addition to AC1-AC3. Certainly, we could improve the reach in this channel, 
however, this is beyond the scope of our analysis that aims to access the impact of the LSP decay on the 
usual SUSY signals. 

For the Oi topology we further veto the presence of isolated leptons with 

p't > 10 GeV and \r^f,\ < 2.5 . (8) 

A charged lepton is considered isolated if the energy deposited in a cone of AR < 0.3 around its direction 
is smaller than 5 GeV. 

The final state 1£ is obtained imposing ([7]) and requiring the presence of only one isolated lepton 
satisfying 

p^T > 20 GeV and \r]e\ < 2.5 . (9) 



9 



Aiming the reduction of the potentially large W background, we further imposed that the transverse 
mass cut 

mT(/^T,PT) > 100 GeV . (10) 

The OS (SS) signal events satisfy ([Tj) and exhibit two isolated leptons with opposite (same) charge. 
The hardest isolated lepton must have 

Pt > 20 GeV and \r]e\ < 2.5 . (11) 

while the second lepton complies with ([5]). 

Trilepton events (3^) must pass the cuts ([7]) and present three isolated leptons with the hardest one 
satisfying (fTTj) and the other two respecting ([8]). Analogously, multilepton events are obtained adding 
one or more isolated lepton to the last topology passing ([5]). 



C. Standard backgrounds 

The main standard model backgrounds for the supersymmetry discovery are 

• QCD processes pp — ^ jjX which have the highest cross section for small values of i?^; 

• tt production that leads to final states WWbb which can lead to many of the final state topologies 
analyzed here. 

• weak gauge boson production in association with jets, denoted by Wj and Zj] 

• production of electroweak gauge boson pairs VV with V = Z or W; 

• Single top production. Here, we did not include the gluon-VF fusion contribution since it is not 
available in PYTHIA. 

We depict in Figure [3] the size and composition of the main SM backgrounds as a function of E^^ for 
SUSY searches. Certainly the SM background for the fully inclusive signal (IN) is the largest one. For 
~ 100 GeV it receives a large contribution from QCD processes leading to light quarks and gluon, 
however, this background decays rapidly with the increase of E^. The main backgrounds for the IN 
topology are tt pair production followed by Wj and Zj productions. At this point, it is important to 
keep in mind that PYTHIA and ISAJET predictions may differ by large factors since they have different 
hypothesis in the fragmentation process, making difficult the comparison of our results with the ones in 



reference 4j| . For this reason and in order to have a consistent calculation of signal and backgrounds we 
have redone all background calculations, instead of using the results of this reference. 

The 0£ final state receives a SM backgrounds similar in size and composition to the IN topology, 
therefore, we do not show the results for this case. We can see from Figure [3] that the bulk of the SM 
background to the 1£ topology originates from tt and Wj production over the entire E^ range. The 
SM background for the OS background is sizeable for E^ < 400 GeV and it is also dominated by the 
ubiquitous tt production as shown in Figure [H The SS background receives its largest contribution from 
tt pair production with one of them coming from the semi-leptonic b decay, however, it is much smaller 
than the previous ones and negligible for E^j. ^ 200 GeV. Furthermore, we verified that the Standard 
Model processes leading to three or more isolated leptons are strongly suppressed, so that the searches 
are basically background free for the E!^ range considered in this paper. 
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FIG. 3: Composition of the inclusive (top left panel), one lepton (top right panel), opposite sign (lower panel) 
backgrounds. The top (red) solid lines stands for the total backgrounds while the blue dashed line corresponds 
to the tt background, the magenta dashed lines corresponds to Wj, and the magenta dotted lines stands for the 
Zj process. The QCD background is represented by the dotted blue line while the dotted black lines stands for 
the single top production and VV is represented by the lower dotted blue line at small E^. 



D. Displaced vertex analysis 

Many of the lightest neutralino decay modes presented in Section can be reconstructed by ATLAS 
and CMS. We concentrated our analysis on the ££, jj, rjj, bb, tt, and t£ modes. We did not try to 
fully identify these decays modes because we are only interested in the possible detection of the displaced 
vertices and not on the further information that can be obtained from them; this analysis will be the 



subject of another work 45[. For instance, in the decays containing r's, we kept only the one- and 
three-prong decays, while we did not study strategies to differentiate the bb from jj modes. In order to 
assure an almost background free search, we required two reconstructed displaced vertices for each event 
that pass the requirements presented below. 

In our analysis we assumed a toy detector based on ATLAS technical proposal [46|. In order to be 
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sure that the LSP decayed away from the primary vertex, we required its decay vertex to be outside an 
ehipsoid 




where the z-axis is along the beam direction. We were conservative on the choice of the size of axis of 
the eUipsoid demanding them to be five times the ATLAS expected errors in each direction - that is, 
Sxy = 20 fim and = 500 //m. We further required that the LSP decay vertex must be inside the pixel 
inner detector, i.e. within a radius of 550 mm and z-axis length of 800 mm. Moreover, to guarantee 
a high efficiency in the reconstruction of the displaced vertices without a full detector simulation, we 
restricted our analysis to reconstructed vertex with pseudo-rapidities \r]\ < 2.5. 

In LSP decays containing r's or b quarks, there might be additional displaced vertices present further 
away from the interaction point. Therefore, to reconstruct the LSP decay vertices we required that all 
visible tracks coming from lightest ncutralino decay cascade must cross a sphere of 10 /im around the 
LSP vertex. 

Relatively long lived SM particles, like r's or B's, can also give rise to displaced vertices as well. 
However, the Standard Model physics background can be eliminated by requiring that the set of tracks 
defining a displaced vertex must have an invariant mass larger than 20 GeV. This way the displaced 
vertex signal passing all the above cuts is essentially background free, except for possible instrument 
backgrounds which are beyond the scope of this analysis. 

At this point we must make sure that the events presenting LSP displaced vertices pass at least of the 
triggers to guarantee that they will be properly recorded. In order to mimic the triggers used by the 
LHC collaborations, we accept events passing at least one of the following requirements: 

• the event has one isolated electron with px > 20 GeV and jryj < 2.5; 

• the event has one isolated muon with pT > 6 GeV and \ri\ < 2.5; 

• the event has two isolated electrons with px > 15 GeV; 

• the event has one jet with px > 100 GeV; 

• the event has missing transversal energy in excess of 100 GeV. 

IV. RESULTS 
A. Standard Analysis 

We estimated the SUSY discovery potential of the LHC in the IN, Oi, U, 3£, M^, OS, and SS 
channels. We analyzed each channel independently, not trying to combine the outcomes of the different 
final state topologies. Certainly combining different topologies enhances the LHC reach for SUSY. For 
the sake of definiteness we defined that a point in the parameter space of our model can be observed 
in a given channel if there is a choice of E"^ for which either the signal leads to 5a departure from the 
background expectation, where this is not vanishing, or 5 events in regions where the SM background 
vanishes. We present our results in the mi/2 ^ mo plane for tan/3 = 10, /i > 0, Aq = —100 GeV, and 
integrated luminosities of 10 and 100 fb~^. 

We depict in Figure H] the LHC discovery potential in the inclusive (left panel) and one lepton (right 
panel) channels assuming R-parity conservation as well as bilinear R-parity violation for an integrated 
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luminosity of 10 fb~^; the results for 100 fb~^ are shown in Figure [5l Current direct search limits coming 
from LEP and Tevatron do not appear in this plot since these constrain lower values of mi/ 2 than the 
ones displayed, highlighting the increased coverage of the LHC. From the left panel of Fig. [Hone can see 
how bilinear R-parity violation reduces the reach in the inclusive channel for a given value of mg with the 
signal depletion growing as mo increases. Basically, the decay of the LSP in BRpV-mSUGRA reduces 
the available missing transverse energy making it harder for the signal to pass the cuts, and consequently, 
stands out from the SM background. This effect is specially important at moderate and high mg where 
the fraction of decays producing neutrinos diminishes; see Fig.[2l Moreover, the shaded area indicates the 
region where the stau is the LSP. As already discussed, this region is excluded if R-parity is conserved 
but not in our BRpV-mSUGRA model. We note that the inclusive signal is rather strong in this area. 
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FIG. 4: LHC discovery reach in the inclusive channel (left panel) and in the one lepton channel (right panel). 
We assumed tan/3 = 10, > 0, Ao = —100 GeV, and an integrated luminosity of 10 fb~^. The points leading 
to visible signal at the LHC are marked by a square for (R-parity-conserving) mSUGRA and by a star for the 
BRpV-mSUGRA case. The (green) shaded area indicates the region where the stau is the LSP. Points already 
excluded by LEP and Tevatron searches lie below the m-1/2 values depicted in this figure. 

The reach in the It channel is also presented in Fig. Hj right panel, for both cases, with and without 
R-parity conservation. Clearly from Fig. 0] one can see that, as expected, the reach in this channel is 
smaller than the one in the inclusive channel with or without R-parity conservation since it diminishes 
as the slepton masses (mo) increase. In this channel, the suppression of the signal is larger than the 
inclusive case when the R-parity violating interactions are present, being visible only at small mo (mi/2) 
in the BRpV scenario. This suppression of the one-lepton BRpV signal is due to the reduction of the 
missing transverse energy at large toq, while at small toq it originates from the production of further 
leptons in the LSP decay as well. Remember that in this channel we require one and only one isolated 
charged lepton. 

The effect of a larger data sample can be seen in Fig. [5] where we assumed an integrated luminosity 
of 100 fb^^. The presence of BRpV reduces the signal reach in the inclusive and one-lepton channels 
when compared to the R-parity conserving case. Moreover, at small and moderate toq the inclusive and 
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FIG. 5: Same as Fig. |4]but for an integrated luminosity of 100 fb ^. 



one-lepton channels lead to similar reach in our BRpV model, however, the inclusive signal has a larger 
reach for mo ^ 1 TeV. 

In Figures [6] and [7] we present the LHC discovery potential in the opposite (OS) and same-sign (SS) di- 
lepton signals (OS (left) and SS (right)) with/without R-parity conservation for integrated luminosities 
of 10 and 100 fb~^ respectively. Clearly, the presence of R-parity violation extends the LHC reach in this 
channel at small toq with respect to the R-parity conserving scenario. On the other hand the coverage 
at large mg is similar in both cases. Again, this extended reach in the BRpV case is due to the presence 
of leptons coming from the LSP decay mediated by light sleptons. In the BRpV scenario the SS channel 
provides the largest reach at moderate and large toq improving the discovery potential while its reach is 
similar to the all inclusive channel for BRpV at small mo. The effect of having a larger data sample, for 
instance 100 fb~^, is shown in Fig. [71 where we can see that the SS signal is more suitable than the fully 
inclusive signal to search for BRpV models at large mo. 

Figs. [S] and [5] display the LHC reach in the three- (left panel) and multilepton (right -panel) channels 
with/without R-parity conservation for integrated luminosities of 10 and 100 fb~^, respectively. The 
discovery potential in these channels is rather limited in the R-parity conserving scenario and it is 
enhanced when we add the BRpV interactions. It is interesting to notice that the poor reach of the M£ 
topology at small mi/2 (^ 200 GeV) for an integrated luminosity of 10 fb~^ is due to the rather hard 
values of used in our analysis. Allowing = 100 GcV leads to a clear at small mi/2. 

For low integrated luminosities (10 fb^^), the 3£ reach is smaller than that of the inclusive channel at 
small and moderate mo, however, it leads to a larger coverage at large mo in BRpV models. This situation 
changes drastically at higher luminosities with the trilepton signal becoming by far the main channel for 
the search of supersymmetry with BRpV. The multilepton channel also becomes very important, being 
the channel with the second largest reach. 

In brief, the decay of the lightest neutralino has an impact of reducing the LHC reach for supersym- 
metry when we take into account only the main all inclusive channel. Notwithstanding, the new BRpV 
interactions lead to a substantial increase in the final states containing isolated charged leptons, specially 
for pairs of same sign leptons and trileptons, compensating the losses in the inclusive channel. If super- 
symmetry is discovered in a set of the topologies discussed above it must be a rather straightforward 
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FIG. 7: Same as Fig. [6] but for an integrated luminosity of 100 fb~^ 



job to establish the presence of R-parity violating interactions. Now, we turn our attention to another 
key characteristic feature of our BRpV-mSUGRA model, namely, the emergence of displaced vertices 
associated to the LSP decay. 



B. Displaced Vertices Analysis 



Following the criteria previously used to define that a signal can be observed in the canonical super- 
symmetry topologies, we have searched for points of the parameter space presenting 5 or more identified 
pairs of displaced vertices for integrated luminosities of 10 and 100 fb~^, assuming that there is no SM 
background for this search. Our results for the LHC discovery reach in the displaced vertex channel are 



15 



^ 1000 
> 



PR conserving trilepton^ 
BRpVtrileptonQ . 

★a ★ 

★u ★ ★ 

ESQ ★ ★ ★ a 

aa a a a a ★ ★ 

aa BaHaa****Ha 

□a aaaaaaaaaaa 

, \ , I , \ , I , 

500 1000 1500 2000 2500 

m„ (GeV) 



^ 1000 

> 



R conserving multileptons^ 

BRpV muitileptonsa 



★ 

★ ★ ★ ★ 

★ ★ ★ ★ ★ 

a* a ★ ★ ★ ★ ★ ★ 
★ 



500 



1000 



1500 



2000 2500 

m„ (GeV) 



FIG. 8: LHC discovery potential in the three lepton channel (left panel) and the multilepton one (right panel) 
for the parameters used in Fig. |4l 



o 

Q 1 ?00 




★ 
a 


★ 








R con 


11 1 1 1 1 . . 

erving trilepto 
BRpV trilepto 


" ■ 

" a 




> 
o 

i 


1?00 












1 ' 


R 


conserving m 
BRpV rr 


ultileptonsB 

ultileptonsQ 


1000 
800 


★ * 

★ i3 

★ Q 


a 
★ 
a 


★ 
★ 
★ 


★ 
★ 
★ 


★ 
★ 


★ 
a 


★ 












1000 

800 




i 

•k-k 


★ 
★ 
★ 


★ 














★ a 


a 


a 


a 


a 


a 


a 


a 


★ 












★ a 


★ 


* 


★ 


★ 








600 


\:a 


a 


a 


a 


a 


a 


a 


a 


a 


a 






600 




*a 


★ 


★ 


★ 


★ 


★ 


★ 






aa 


a 


a 






a 






a 












a a 


a 


a 


a 


a 


★ 


★ ★ 


★ 


400 


aa 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 




400 




aa 


a 


a 


a 


a 


a 


★ ★ 






aa 


a 

, , , 


a 


a 

, , 


a 


a 


a 

, 


a 


a 




a 


a 








aa 


a 


a 


a 


a 

1 , 


a 


a ★ 


a a a 



250 500 750 1000 1250 1500 1750 2000 2250 250 500 750 1000 1250 1500 1750 2000 2250 



m„ (GeV) m„ (GeV) 

FIG. 9: Same as Fig. [8] but for an integrated luminosity of 100 fb~^. 

also shown in the toq ® rai/2 plane for tan/3 = 10, /i > and Aq ~ —100 GeV. 

In Figure [To] we present the displaced vertex reach. As one can see from this figure, the LHC will be 
able to look for the displaced vertex signal up to mi/2 ~800 (1000) GeV {rn^o ^ 340 (430) GeV) for a 
large range of toq values and an integrated luminosity of 10 (100) fb~^. Notice that the reach in this 
channel is rather independent of mg as expected from Fig. [TJ Moreover, this signal for BRpV-mSUGRA 
disappears in the region where the stau is the LSP due to its rather short lifetime. It is interesting to 
compare these results with ones obtained in the canonical search channels. The displaced vertex signal 
provides the largest reach at moderate and large toq, being more important that the IN, and 
topologies. Notwithstanding its good reach at small toq, it is not the main search channel in this range. 
Certainly, the combination of the different final state topologies will help us to pinpoint the properties of 
the model. Furthermore, we should not forget that the discovery of high invariant mass displaced vertices 
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allow us to study directly the LSP such as its mass and branching ratios. Moreover, we verified that 
displaced vertices possess a larger reach than the standard search channels not only at large mo, but also 
at large mi/2- We expect that for larger values of tan/3 the reach in the displaced vertex analysis will be 
even better 

Since the LHCb experiment has excellent vertex capabilities, it is natural to conclude that the LHCb 
can also play an important part in the displaced vertex search for supersymmetry, at least in the first 
stages of LHC operation, i. e. when the accumulated luminosity is rather small [45j . 



V. CONCLUSIONS 

We have considered bilinear R~parity violating supergravity, the simplest effective model for the super- 
symmetric origin of neutrino masses, reproducing the current neutrino oscillation data. We have given a 
detailed analysis of the impact of R-parity violation on the standard supersymmetry searches. We have 
examined how the LSP decays via bilinear R-parity violating interactions degrade the average expected 
missing momentum in the reactions and shown how this diminishes the reach in the main canonical SUSY 
search channel that is the all inclusive topology. However we have seen how this can be overcompensated 
by the coverage obtained through the search for final states containing isolated charged leptons, specially 
the trilepton signal. Moreover, we can also search for displaced LSP decays which can further increase the 
LHC reach for SUSY. We present the resulting reach in supergravity parameter space for two reference 
luminosities of 10 & 100 fb"^ 

In the first stages of LHC operation, when the accumulated luminosity is rather small, we stress that 
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a displaced-vertex analysis by the LHCb collaboration can play an important role in probing bilinear R- 
parity violating supergravity. With enough luminosity to be collected in subsequent stages of operation, 
the challenge will be to analyze the detailed nature of the LSP decays as a test of the neutrino mixing 
angles measured in oscillation experiments. We plan to turn to these issues elsewhere. 
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